2 



Strain induced magnetic domain evolution and spin re-orientation 
transition in epitaxial manganite films 

Gyanendra Singh J*| P. K. RoutJ*| and Rajni Porwal 

Condensed Matter - Low Dimensional Systems Laboratory, Department of Physics, 
Indian Institute of Technology Kanpur, Kanpur - 208016, India 

R. C. Budhanfl 

Condensed Matter - Low Dimensional Systems Laboratory, Department of Physics, 
Indian Institute of Technology Kanpur, Kanpur - 208016, India and 
National Physical Laboratory, New Delhi- 110012 India 

Abstract 

The evolution of magnetic domain structure in epitaxial Lao.62sCao.375Mn03 films on (001) 
NdGa03 is monitored as a function of temperature and magnetic field using Magnetic Force Mi- 
croscopy. We see two distinct regions of magnetic orientational order; one in-plane displaying 
contrast-less image and the other tilted away from the film plane forming a distinct stripe pattern. 
A strong domain splitting is observed at the boundary of two regions, which is resilient to reorien- 
tation with temperature and magnetic field. We propose a model magnetic free energy functional 
to explain the mechanism of domain splitting seen in manganite films. 
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Magnetic domains (MDs) in ferromagnetic films arise from the requirement of mini- 
mization of the total magnetic free energy consisting of magnetic interactions of both local 
and non-local nature. 1 The epitaxial films of Lai_ x Sr x Mn03 (LSMO) and Lai_ x Ca x Mn03 
(LCMO) provide suitable platforms to see rich magnetic textures because of the coupling 
between the charge, spin, orbital and lattice degree of freedom, 2 which connect to the various 
components of the magnetic free energy functional, and are affected significantly by epitaxial 
strain. The latter is a powerful tuning parameter as these oxides can be grown epitaxially on 
a large number of single crystal substrates which impart varying degree of strain, depending 
on the lattice mismatch. Due to the magnetoelastic coupling these elastic strains can induce 
magnetic anisotropy in the film, whose magnitude depends upon the magnetostriction con- 
stants and the amount of stress in the film. Typically the value of stress induced anisotropy 
lies within 10 4 -10 5 Jm -3 for thin films of various manganites.- - - Moreover the intrinsic mag- 
netocrystalline anisotropy is in the order of ~ 10 3 Jm -3 , which is quite small as compared 
to the stress induced anisotropy.- Thus the strain plays an important role in determining 
the magnetic properties of the epitaxial films. The compounds LCMO and LSMO are pseu- 
docubic perovskites with lattice parameters in the range of 0.386-0.389 nm in the unstrained 
form. The commonly used substrates for epitaxial growth of these oxides are LaA103 (LAO), 
SrTi03 (STO) and NdGa03 (NGO). While LAO provides an isotropic in-plane compressive 
strain, which results in out-of-plane easy axis and maze-like domains, the films on STO have 
the easy axis lying in-plane due to in-plane tensile strain and thereby form planar domains .- 
The orthorhombicity of NGO makes the compressive strain anisotropic in the film plane 
and thus promoting a preferential direction for MD formation. However the imaging of MD 
structure of the films on NGO has given contrasting results, with observation of both in- 
plane and out-of-plane orientation of magnetization.- - — The reasons for this non-uniqueness 
of the magnetic texture have not yet been established. 

Here we report a careful study of the evolution of MD structure in La .625Ca .375MnO3 
films as a function of temperature (T) and in-plane magnetic field (H\\) in ultra high vacuum 
of ~ 10 -10 Torr using Magnetic Force Microscopy (MFM) (Scanning Probe Microscope, 
Model: UHV 3500, RHK Technology). We first show the manifestations of the magnetic 
ordering in the magnetization M(T) and resistance R(T) data of a 200 nm thick LCMO film 
grown on (001) NGO at 800°C by pulsed laser ablation of a target of La .625Ca .375MnO3. 
The details of the thin film growth are described in several earlier publications. 13 The Tc 
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as estimated from the R(T) and M(T) data is ~ 260 K [Inset of Fig. 1(a)]. Figure 1(a) 
shows M{H) loops measured with H\\ along [100] direction, as well as in out-of-plane field 
H± configuration. From the shape of the loops, it is evident that M has both in-plane and 
out-of-plane components, with in-plane remanence being ^30% of the saturation value as 
deduced from H\\ loop. 

The MFM image of the film shows a regular stripe domain pattern with the bright 
and dark contrast due to a quasi-periodic orientation of the component of M in to and 
out of the film plane [Fig. 1(b)]. The period of the stripe domains (L) is 300 nm 
as determined by Fast Fourier Transformation (FFT) shown in the inset of Fig. 1(d) of 
the MFM image. Although the homogeneous stripe domain pattern observed here has been 
reported previously, 6 the striking discovery of this study is the observation of the microscopic 
patches of two magnetically distinct phases separated by a thin boundary marked by line B 
in the MFM image shown in Fig. 1(c). A remarkable feature of Fig. 1(c) is the branching 
of the stripes as the boundary of the CL region is approached. This splitting of the domains 
has been quantified by taking line scans perpendicular to the length of the stripes at various 
points across the boundary from line A to B and the variation of L with distance is shown 
in Fig. 1(d). While the stripes multiply, the MFM signal intensity drops on approaching the 
boundary. One may suspect that the CL phase is confined to the thinner sections of the film 
where the dominant effect of dipolar energy will try to make M in-plane. To address this 
possibility, we have done MFM imaging after Ar + ion milling of the film down to 100 nm 
but the stripes persist even at this thickness. Clearly, the effect is not due to a variation in 
film thickness. The simplest explanation for these observations can be based on the spatial 
variation of strain in the film. While a previous report on LCMO grown on (001) NGO 
reveals a critical thickness of ~ 500 nm for fully relaxed bulk like state,— it is possible that 
at intermediate thickness of ~ 200 nm of our film both strained and relaxed regions coexist 
due to partial strain relaxation. Although the branching of domains is expected near the 
edge of the sample, all CL regions observed in the film lie well within the sample boundaries. 
We have taken a number of scans at different places of the film and observed only <~ 5% of 
such CL regions. 

A more informative approach to understand MDs would be to see the variation of the 
angle 9 between M and film plane [Fig. 1(e)], which has been extracted as 6{x ) y) = 
sm- 1 [MFM(x,y)/MFM max } ) where MFM(x,y) is the intensity of the MFM signal and 
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MFM max is its maximum value corresponding to 9 = ±90°. 15 We see two separated regions; 
one with 9 ~ 0° (in-plane M) and the striped region with non zero #'s. Surprisingly, near 
the phase boundary we see a region with maximum canting angle 9 m (dark red and dark 
blue), which is much higher than 9 m « 65° seen in the regions away from the boundary. 
Moreover this 9 m « 65° is distinct in the histogram extracted from the 9 image [Fig. 1(f)]. 
The sinusoidal oscillation of M z with a particular 9 m would give an equal number of counts 
to the 9 values ranging from +9 m to —9 m . 20 This clearly means a symmetric histogram 
profile about 0°. Thus a more common value of the 9 m implies discontinuities at +9 m and 
— 9 m in histogram. We indeed observe discontinuities in 63-67° range on both branches of 
histogram. 

In order to address the micromagnetic domain evolution at the boundary region in the 
vicinity of Tq : the MFM images are taken at various temperatures. At T < 240 K, both 
the phases with domain branching at their boundary are observed (See Fig. 2). As the 
temperature increases to 254 K, the branching disappears, but the stripes remain, though 
with attenuated intensity. On further increasing the temperature to 261 K, the stripe do- 
mains disappear completely in the interior regions; while some magnetic contrast persists at 
the boundary till 264 K. It is clear that the magnetic order parameter at the boundary is 
much more resilient to change, suggesting some kind of a pinning mechanism in action. The 
possible origin of such magnetic inhomogeneity can be found in Ref. [13], where mesoscale 
regions with different Tq are observed near the artificial grain boundaries in LSMO films 
grown on bicrystal STO substrates. They have attributed such effect to the spatial variation 
in the strain, which we believe to be the reason for non-trivial domain pattern seen in our 
film. A similar resilience of boundary to reorient is observed for in-plane magnetic fields up 
to « 100 niT. 

All these observations can be explained by a simple model in which the film lies on the 
x?/-plane and M has both y and z components in addition to M z oscillating as sin(7rx/L) 
as sketched in Fig. 3(a).— ^ The free energy density E of such a domain pattern can be 
expressed as: 



Here the first term is the exchange interaction energy with A = 1.7x10 12 Jm 1 . 19 The 



E = 




(1 - cos 9 m ) + - {Ki cos 2 9 m + K 2 cos 4 9 m ) 



^ 0m _2joHMs I{0m) 



(i) 
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second term is the uniaxial magnetic anisotropy energy taken up to the fourth order and 
the next two terms are the short and long range part of dipolar energy with Q = |/x M| = 
1.01 x 10 5 Jm -3 . The last term is the Zeeman energy for an in-plane magnetic field H applied 



along y-axis with I(8 m ) = Jq /2 dtpyjl- sin 2 9 m sin 2 Lp « 1 + (f - l) (cos# m ) 3/2 . 

Minimization of Eq. (1) relative to L and 9 m yields an analytical expression for domain 
period: 

L = (2) 
f2oo(l + cos# m ) 

where cos 9 m satisfies the relation 

ft 2 a 2 ( VL 2 a 2 \ 

3 ~32A C ° s2 9m + ( Kl + 2X2 C ° s2 ° m ~ Q + ~16A J C ° S ° m 

-^(f-l)(cos^-^0 (3) 

The stripe domains present in our film can be explained by the canted state of M in- 
troduced solely by a non-zero K2 20 The least square fit shown in Fig. 1(a) of the in-plane 
magnetic hysteresis data yields the anisotropy values of Ki/Q ~ 0.95 (or K\ « 9.60 xlO 4 
Jm -3 ) and K 2 /Q ~ 0.12 (or K 2 ~ 1.21 xlO 4 Jm -3 ). 15 In comparison the value of mag- 
netocrystalline anisotropy is at least an order of magnitude smaller. The stress induced 
anisotropy for biaxial strains can be expressed as K s = 3Aa/2, where A is the magnetostric- 
tion constant and a is the stress. 4 The stress can be estimated from the product of the Youngs 
modulus (Y) and the strain (e). The typical values of e ~ 10 -2 and Y ~ 10 n -10 12 Nm" 2 
result in stress values of ~ 10 9 -10 10 Nm" 2 . 5 Using literature value of the magnetostriction 
constant A ~ 10 -4 , 21 one obtains an anisotropy of K s ~ 10 5 Jm -3 , which is quite close to 
our value of K\. This indicates that the anisotropy present in the system is predominantly 
due to the elastic strains. With L ~ 300 nm as calculated before, we have 9 m « 65°, which 
is in agreement with the value determined from 9 image. 

Equation (2) reveals that as 9 m increases from 0° to 90°, L becomes twice as large. 
Thus, as we move from stripe domain region to CL region where 9 m = 0°, one can expect 
a subdivision of the stripes. A simple way to explain the mechanism responsible for the 
change in 9 m will be a variation of magnetic anisotropy due to local elastic strains. The 
minimization of total anisotropy energy along the lines of Ref. [18] under the appropriate 
strain conditions (e xx « e yy ^ and e xy = 0) shows a direct proportionality of X's with 
the in-plane strain, which means higher the compressive strain higher will be the anisotropy 



and vice versa. All the possible magnetic domain configurations are summarized in Fig. 
3(b) for K1/K2 = 7.9 in our film. We can clearly see that Ki/Q = 0.80 is the critical 
anisotropy value separating planar domains from stripe ones. For K\/VL < 0.80 (region- 
I), M is completely in-plane whereas the stripe domains with canted M can be observed 
for Ki/Vt > 0.80 (region-II & III). The maximum canting angle 9 m gradually increases in 
region-II and reaches a value of 90° in region-Ill. A slight decrease or increase in strain from 
the present value (Ki/Q = 0.95) can move the system towards region-I or III respectively 
and thus a transition between planar and stripe domain pattern as observed here can be 
expected. Furthermore the contradicting domain structures reported previously in Ref. [3-9] 
can be due to the strain present in the films depending on the growth conditions and the 
film thickness. Similar to the nature of L, we see a minor increase in K\/VL followed by a 
gradual drop whereas the corresponding # m 's decrease rapidly to 0° from a value of 65° while 
reaching a maximum of 85° at the middle [See Fig. 3(c)]. Although phase separated regions 
resulting from the local variation in strain have been proposed before,— we have given a 
direct and visual evidence of different M states with same magnetic ordering temperature 
in LCMO film. 

In summary, MFM has been used to establish the two distinct orientations of M in 
La .625Ca .375MnO3 epitaxial films. A stripe domain pattern caused by a tilted orientation 
of M with respect to the film plane shows a distinct subdivision in a manganite film. A 
model calculation predicts the fragmentation of the stripes at the boundary due to the local 
variation of strain. 

The authors thank P. C. Joshi for technical help. This research has been supported by 
grants from DIT. G.S. and P.K.R. acknowledge the financial support from CSIR, India. 
R.C.B. acknowledges the J. C. Bose fellowship of DST. 



* G.S. and P.K.R. contributed equally to this work. 
^ rcb@iitk.ac.in, rcb@nplindia.org 

1 A. Hubert and R. Schaffer, Magnetic Domains: The Analysis of Magnetic Microstructures 
(Springer, New York, 1998), p. 99. 

6 



2 Y. Tokura, Rep. Prog. Phys. 69, 797 (2006). 

3 S. E. Lofland, S. M. Bhagat, H. L. Ju, G. C. Xiong, T. Venkatesan, R. L. Greene, and S. Tyagi, 
J. Appl. Phys. 79, 5166 (1996). 

4 J. O'Donnell, M. S. Rzchowski, J. N. Eckstein, and I. Bozovic, Appl. Phys. Lett. 72, 1775 
(1998). 

5 L. M. Berndt, V. Balbarin, and Y. Suzuki, Appl. Phys. Lett. 77, 2903 (2000). 

6 J. Dho, Y. N. Kim, Y. S. Hwang, J. C. Kim, and N. H. Hur, Appl. Phys. Lett. 82, 1434 (2003). 

7 R. Desfeux, S. Bailleul, A. Da Costa, W. Prellier, and A. M. Haghiri-Gosnet, Appl. Phys. Lett. 
78, 3681 (2001). 

8 A. Biswas, M. Rajeswari, R. C. Srivastava, T. Venkatesan, R. L. Greene, Q. Lu, A. L. de 
Lozanne, and A. J. Millis, Phys. Rev. B 63, 184424 (2001). 

9 C. Israel, L. Granja, T. M. Chuang, L. E. Hueso, D. Snchez, J. L. Prieto, P. Levy, A. de Lozanne, 
and N. D. Mathur, Phys. Rev. B 78, 054409 (2008); D. Sanchez, L. E. Hueso, L. Granja, P. 
Levy, and N. D. Mathur, Appl. Phys. Lett. 89, 142509 (2006). 

S Piano, A. D. Santis, F. Bobba, F. Giubileo, M. Longobardi, A. D. Bartolomeo, M. Polichetti, 
A. Scarfato, D. Zola, A. Vecchione, and A. M. Cucolo, J. Phys. Condens. Matter 21, 254205 
(2009). 

1 S. J. Lloyd, N. D. Mathur, J. C. Loudon, and P. A. Midgley, Phys. Rev. B 64, 172407 (2001). 

2 V. K. Vlasko-Vlasov, Y. Lin, U. Welp, G. W. Crabtree, D. J. Miller, and V. I. Nikitenko, J. 
Appl. Phys. 87, 5828 (2000). 

3 P. Padhan and R. C. Budhani, Phys. Rev. B 67, 024414 (2003); P. K. Muduli, G. Singh, R. 
Sharma, and R. C. Budhani, J. Appl. Phys. 105, 113910 (2009). 

4 S. Seiro, E. Roller, Y. Fasano, and 0. Fischer, Appl. Phys. Lett. 91, 091913 (2007). 

5 See supplementary material at [URL will be inserted by AIP] for the theoretical analysis of 
MFM images and the least square fitting of hysteresis data. 

6 Y.-A. Soh, G. Aeppli, N. D. Mathur, and M. G. Blamire, Phys. Rev. B 63, 020402(R) (2000). 

7 Y. Yafet and E. M. Gyorgy, Phys. Rev. B 38, 9145 (1988). 

8 Y. Z. Wu, C. Won, A. Scholl, A. Doran, H. W. Zhao, X. F. Jin, and Z. Q. Qiu, Phys. Rev. Lett. 
93, 117205 (2004). 

9 J. W. Lynn, R. W. Erwin, J. A. Borchers, Q. Huang, A. Santoro, J. L. Peng, and Z. Y. Li, 
Phys. Rev. Lett. 76, 4046 (1996). 

7 



In the absence of the fourth order term {K2 =0), using the definition of the effective anisotropy 
= \ J (M_l — M\\)dH, we find K^^ /CI « —0.34, which suggests a homogeneous in-plane 
magnetization description and thus a planar magnetic domain structure. 

M. R. Ibarra, P. A. Algarabel, C. Marquina, J. Blasco, and J. Garcia, Phys. Rev. Lett. 75, 3541 
(1995). 

A. E. Clark, Ferromagnetic Materials edited by E. P. Wohlfarth (North-Holland, Amsterdam, 
1980), p.531. 

A. I. Tovstolytkin, V. V. Dzyublyuk, D. I. Podyalovskii, X. Moya, C. Israel, D. Snchez, M. E. 
Vickers, and N. D. Mathur, Phys. Rev. B 83, 184404 (2011). 



8 



1.0 
0.5 

CO 

0.0 
-0.5 



-1.0 



(a) 

[001] 

-h*— H|| [100] 






^10 
|o.5 
^0.0, 


100 30C 





-100 -50 



H (mT) 



50 



100 





0.0 0.5 1.0 1.5 2.0 

Distance (//im) 




-90 -60 -30 30 60 90 

6 (Degree) 



FIG. 1. (a) Hysteresis loops measured with field (H) along [100] and [001] directions. The solid 
line shows the least square fit to in-plane hysteresis loop. 15 The inset shows the R vs. T along with 
zero field cooled (open circle) and field cooled (solid circle) M vs. T measured in H = 20 mT. 
(b & c) The MFM images of LCMO film taken at two different places at 120 K in the absence of 
magnetic field after zero field cooling. Here the scan area is 4x4 /im 2 . (d) The domain period as 
a function of length measured from line A to B as shown in (c) . The inset shows the FFT image 
of the panel (b). (e) The 9 image extracted from MFM image (c). (f) The histogram of 9 image 
showing most probable 9 m values (blue regions). 
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FIG. 2. (a-f) The MFM images at several temperatures near Tq in zero magnetic field. All images 
captured have the scan area of 4x4 /im 2 taken at the same place in continuous heating mode while 
the nanometer scale offset due to the thermal drift has been controlled by x-y offset. 
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FIG. 3. (a) Schematic image the stripe domain pattern with M lying completely in ?/z-plane with 
M z (x) ~ sin(7rx/L). The right panel show the MFM signals as we move across the stripe length, 
which fits quite well with a sinusoidal function (red line), (b) The 6 m and L as a function of K±/Q 
with K1/K2 = 7.9. The projections of M on yz-plane are shown for different regions, (c) The 
variation of L, Ki/ft and 9 m as we move from line A to B shown in Fig. 1(c). 
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